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ABSTRACT 

The diffuse cm-wave IR-correfated signaf, tfie "anomafous" CMB foreground, is 
thought to arise in the dust in cirrus clouds. We present Cosmic Background Im- 
ager (CBI) cm-wave data of two translucent clouds, ^ Oph and LDN 1780 with the 
aim of characterising the anomalous emission in the translucent cloud environment. 

In C, Oph, the measured brightness at 31 GHz is 2.4cr higher than an extrapolation 
from 5 GHz measurements assuming a free-free spectrum on 8 arcmin scales. The SED 
of this cloud on angular scales of 1° is dominated by free-free emission in the cm-range. 
In LDN 1780 we detected a 3 cr excess in the SED on angular scales of 1° that can be 
fitted using a spinning dust model. In this cloud, there is a spatial correlation between 
the CBI data and IR images, which trace dust. The correlation is better with near-IR 
templates {IRAS 12 and 25 /im) than with IRAS 100 /^m, which suggests a very small 
grain origin for the emission at 31 GHz. 

We calculated the 31 GHz cmissivitics in both clouds. They are similar and have 
intermediate values between that of cirrus clouds and dark clouds. Nevertheless, we 
found an indication of an inverse relationship between emissivity and column density, 
which further supports the VSGs origin for the cm-emission since the proportion of 
big relative to small grains is smaller in diffuse clouds. 

Key words: ISM: individual(LDN 1780, C, Ophiuchi) - radiation mechanism: general 
- radio continuum: ISM 



1 INTRODUCTION 

Amongst the challenges involved in the study of Cosmic Mi- 
crowave Background (CMB) anisotropies is the subtraction 
of Galactic foreground emission. The detailed study of these 
foregrounds led to the discovery of a new radio-continuum 
mechanism in the diffuse interstellar medium (ISM). It was 
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first detected by the Cosmic Background Explorer as a dif- 
fuse all-sky dust-correlated si gnal with a flat sp ectral in- 
dex between 31 and 53 GHz ('Kogut et al."l996a||3). With 
ground based observations, Lcitch ct al. ( 1997|) observed 
high Galactic latitude clouds and detected an excess at 
14.5 GHz that, because of a lack of Ha emission, could 
not be accounted for by free-free with typical temperatures 
(T - lO'^K). 

Different emission mechanisms have been pro- 
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Figure 1. {u,v) coverage of the CBI (left) and CBI2 (right) in the configurations used for the observations of Oph on the left and 
LDN 1780 on the right. Axes are in units of wavelength (A). 



posed for the anomalous 
dust IP raine fc LazarianI 
Draine fc Lazarian. ,1999, ). 
Leitch et all Il997l) and 
Bennett et al] |2003| ). To 



emissio n, such as spinning 
Il998al lbl). magnetic dust 
hot (T~10'' K) free-free 
hard synchrotron radiation 
date, the evidence favours 



the spinning dust gra in mode l jFinkbeiner et al] 



de Olive ira-Costa et al. I |2004|; I Watson et al.l 



2004 



2005 



"Casassus et al.1 |2006| . l2008l : ICastellanos et al.l |2010| ) in 
which very small grains (VSGs) with a non-zero dipole 
moment rotating at GHz frequencies emit cm-wave ra- 
diation. The emission has its peak at ~20-40 GHz and 
the model predicts that it is dominated by the smallest 
grains, possibly po ly cyclic aromatic hydrocarbons (PAHs). 
lYsard et al.l (|20ld ) found a correlation across the whole sky 
using WMAP data between 23 GHz maps and IRAS 12 /.im, 
which supports the VSG origin for the cm-wave emission. 

The radio-IR correlation suggests that the "cirrus" 
clouds are responsible for the anomalous emission at high 
Galactic latitudes. Cirrus is the large-scale filamentary 
structure detected by IRAS ||Low et al.l Il984l ). It is seen 
predominantly at 60 and 100 fim in the IRAS bands and 
the origin of the IR-radiation is generally ascribed to dust 
continuum emis sion with a contribution from atomic lines 
(e.g. Ol 63 um. IStar^ll990l ). On shorter wavelengths (e.g. 
12 fim), the emission is from polycyclic aromatic hydrocar- 
bons (PAHs) subject to thermal fluctuations. Cirrus clouds 
span a wide range of physica l parameters ; most arc a tomic 
and some are molecul a r (e.g. Miville-Deschenes et al.l l2002l : 
ISnow fc McCaill2006l ). IStark et al.l (|l994h . using a combi- 
nation of absorption and emission line measurements, found 
Hi gas temperatures in the range between 20 and > 350 K. 
The low temperatures correspond to the coldest clumps in 
the clouds that represent the cores of the much more widely 
distributed and hotter Hi gas. Cirrus clouds are pervaded by 
the interstellar radiation field (IRF). They have column den- 
sities A^(Hl)~ 1 — lOx 10^*^ cm~^, which correspond to visual 



extinctions of Ay ^ 1 mag assuming a typical dust-to-gas 
ratio of 100. It is commonly found that these clouds are not 
gravitationally boun d and that their kin ematics are domi- 
nated by turbulence (|Magnani et al.iri985l ). Cirrus clouds are 
difficult to characterise because of their low column densi- 
ties, and are thus not the best place to study the anomalous 



emission. 



Local known clouds offer an opportunity to char- 
acterise the cm-emission. A number of investigations 
have been made of on these c louds: the Perseus molec- 
ular cloud l| Watson et al.l |2005|). o Oph JCasassus et al.l 
l2008l). RCW 175 jPickinson et ahl |2009|). LDNllll 
ll Ami Consortium et al.l |2009|). M T 'S (ICastellanos et al.1 
l2010l ). among others. iFinkbeineJ (|2004h using WMAP 
(23-94 GHz) and Green Bank (5,8,10 GHz) data fitted 
a spinning dust model to the e mission from the cloud 
LDN1622. ICasassus et al.l l|2006l ). with interferometric 
observations at 31 GHz of the same cloud, detected bright 
cm-emission where no emission from any known mechanism 
was expected. They performed a morphological analysis 
and found a better cross-correlation betwee n 31 GHz and 
IRAS 12 fim than with IRAS 100 /im. IWatson et al] 
(12005 ) with the COSMOSOMAS experiment discovered 
strong emission in the frequency range 11-17 GHz in the 
Perseus Molecular Cloud (G159.6-18.5). In this cloud, 
the spectral energy distribution (SED) is well-fitted by a 
spinning dust model. A disadvantage of these studies, in the 
aim of characterising anomalous emission as a foreground 
to the CMB, is that the clouds already studied potentially 
constitute a very different phase of the ISM than the cirrus, 
so inferences drawn from them may not be applicable to 
understanding the anomalous emission from the cirrus seen 
at high Galactic latitudes. 

In this paper we present cm- wave continuum data of the 
translucent clouds Oph and LDN 1780 acquired with the 
Cosmic Background Imager (CBI). Translucent clouds are 
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Figure 2. CLEANed 31 GHz images of C Oph (left) and LDN 1780 (right). Contours are 10, 30, 50, 70, 90% of the peak brightness, 
which is 9 X 10"^ MJy sr'^ in C Oph and 0.016 MJy sr^^ in LDN 1780. The FWHM primary beams of CBI and CBI2 are shown as 
dashed lines (in LDN 1780, there are two pointings). These images have not been beam-corrected. The point source NVSS 153909-065843 
in LDN 1780 is marked with an arrow. The synthesised beam is shown at the bottom right corner of each frame. 



interstellar clouds with some protection from the radiation 
field in that their extinction is in the range Av ~ 1 — 4 mag 
jSnow fc McCaljlioOS ). They can be understood as photo- 
dissociation-regions (PDR). Across translucent clouds, car- 
bon undergoes a transition from singly-ionised into neutral 
atomic or molecular (CO) form. In translucent clouds, physi- 
cal properties such as density and temperature, and environ- 
mental conditions, such as exposure to the IRF, are interme- 
diate between those of dense clouds and those of transparent 
(cirrus) clouds. By bridging the gap in physical conditions, 
translucent clouds are test beds for the extrapolation of the 
radio/IR relative emissivities seen in dense clouds. 

The Lynds Dark Nebula (LDN) 1780 is a high Galac- 
tic latitude (/ = 359.0°, b = 36.7°) tra n slucent region at 
a distance of llOilO pc (|Francd [l989l ) . iRidderstad et al.1 
found that the spatial distribution of the mid-IR 
emission differs significantly from the emission in the far- 
IR. Also, they show using IR colour ratios that there is an 
overabundance of PAHs and V SGs with respect to the sola r 
neighbourhood (as tabulated in lBoulanger &: PeraultlfigSSl ). 
although their result can be explained by an IRF that is 
overabundant in U V photons compared to the standard IRF 
(|Witt et al.ll2010l ). Using an opt ical-depth map constructed 
from ISO 200 ^m emission, [Riddcrstad ct al. (2006) found 
a mass of ~ 18 Mq and reported no young stellar objects 
based on the absence of colour excess in point sources. Be- 
cause of the morphological differences in the IR, this cloud 
is an interesting target to make a morphological comparison 
with the radio data, in order to determine the origin of the 
anomalous emission. The free-free emission from this cloud 
is very low, which is favourable for a study of the radio-IR 
correlated emission. 

The cloud coincident with ( Oph is a prototypical and 



well-studied translucent clou d. C Oph itself is an 9Vb star 
at a distance of 140±16 pc (fPerrvman et al.lll997l ). In this 
line of sight, the total H-nucleus column density is fairly 
well determined, iV(H) ~ 1.4y_l^cm-^, with 56% of 
the nuclei in molecular form (|Mortonlll975l ). Although the 
observations reveal several interstellar components at dif- 
ferent heliocentric velocities, the one at Vq — —14.4 km 
s~^ contains most of the material and is referred to as 
the Oph cloud. There have be en numerous efforts to 
build chemical models of this cloud ("Black fc Dalgarnolll977l : 
Ivan D ishocck fc Black 198J)i^iala et al. 198^). In the pio- 
neering work of lBlack fc Dalgarnol (Il977l ). a two shell model 
is proposed to fit the observational data: a cold and denser 
core surrounded by a diffuse envelope. The density and tem- 
perature in these models are in the range n — 250 — 2500 
cm"^, T = 20 - 100 K for the core and n = 200 - 500 
cm""^, T = 100 — 200 K for the envelope. These conditions 
approach those of the cirrus clo ud cores whose densities lie 
in the same range (|Turnej|l99l ). 



The rest of the paper is organised as follows: in § 2 
we describe the data acquisition, image reconstruction and 
list the auxiliary data used for comparison with the CBI 
data. § 3 contains a spectral and morphological analysis of 
both clouds. In § 4 there is a comparison of the 31 GHz 
emissivity of dark, translucent and cirrus clouds. § 5 presents 
our conclusions. 
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2 DATA 



2.1 31 GHz Observations 

The observations at 31 GHz we re carried out with the Cos- 
mic Background Imager (CBIl. IPadin et al]|2002l ). a 13 el- 
ement interferometer located at an altitude of 5000 m in 
the Chajnantor plateau in Chile. Each antenna is 0.9 m in 
diameter and the whole array is mounted on a tracking plat- 
form, which rotates in parallactic angle to provide uniform 
uv coverage. The primary beam has full-width half maxi- 
mum (FWHM) of 45.2 arcmin at 31 GHz and the angular 
resolution is ^^8 arcmin. The receivers operate in ten fre- 
quency channels from 26 to 36 GHz. Each receiver measures 
either left (L) or right (R) circular polarisation. The inter- 
ferometer was upgraded during 2006-2007 with 1.4 m dishes 
(CBI2) to increase temperature sensitivity (Taylor et al. in 
prep.). The primary beam is 28.2 arcmin FWHM at 31 GHz 
and the angular resolution was increased to ^^4 arcmin. 

Table [T] summarises the observations. Both clouds were 
observed in total intensity mode (all receivers measuring L 
only). C, Oph was observed using CBIl in a single position 
while LDN 1780 was observed in two different pointings: 
L1780E and L1780W by CBI2. The configurations of the 
CBIl and CBI2 interferometers result in the (u, v) coverage 
shown in Fig. [T] 

We reduced the data using th e same routin es as 

those used for CMB data analysis jPearson et al.l l2003l : 
iReadhead et all l2004al lbl). Integrations of 8-min on source 
were accompanied by a trail field, with an offset of 8 arcmin 
in RA, observed at the same hour angle for the subtraction 
of local correlated emission (e.g. ground spillover). Flux cal- 
ibration is tied to Jupiter (with a brightness temperature of 
146.6 ±0.75 K. lHill et al1l2009l ). 



Table 1. Summary of observations. The coordinates are J2000 
and the time quoted is the observation time on-souree. 



Source Date R.A. dccl. Time 

Oph 8/7/04 16'^37™9=.5 -10°34'01" 8000 s 

2L1780E 13, 28/4/07 15'M0™30= -07° 14' 18" 12000 s 

2L1780W 17, 18/4/07 15'^39™40= -07°11'40" 8000 s 



^: observed with CBIl 
observed with CBI2 



2.1.1 Image reconstruction 

The reduced and calib rated visibilitie s were imaged with 
the CLE AN algorithm (|Hogbomlll974h using the DIFMAP 
package l|Shepherdlll997f ). We chose natural weights in or- 
der to obtain a deeper restored image. The theoretical noise 
level (using natural weights, as expected from the visibility 
weights evaluated from the scatter of individual frames) is 
4.9 mjy beam"^ for C Oph. Fig. [2] shows the CLEAN im- 
age of Oph. This image was then corrected by the pri- 
mary beam response of CBIl at 31 GHz (45.2 arcmin). 
For LDN 1780, we obtained CLEAN images of the two 
fields. The estimated noise level is 4.3 mJy beam~^ for 
L1780E and 3.0 mJy beam"^ for L17 80W. In LDN 1780W 
a point source, NVSS 153909-065843 (|Condon et al.lll998l '). 
lies close to the NNW sector of the half-maximum con- 
tour of the CBI2 primary beam in Fig. (2] This source al- 
lowed us to fix the astrometry of the data which corre- 
sponded to an offset of 30" and 56" in R.A. and Dec. These 
shifts are consistent with the root mean square (rms) point- 
ing accuracy of ~ 0.5 arcmin. The restored CBI2 fields 
of LDN 1780 were then combined into the weighted mo- 
saic shown on Fig. 2, according to the following formula: 
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Table 2. Auxiliary data used. 



Survey /Telescope 


X/i/ 


6*1 


Reference 


SHASSA 


656.3 nm 


0.8 arcmin 


Gaustad et al. f2001) 


Spitzer^ 


8 


2" 


Fazio et al. (2004) 


IRAS /mis 


12, 25, 60 & 100 fj.m 


3.8 - 4.3 arcmin 


Miville-Deschenes & Lagaclie (2005) 


ISO 


100 & 200 fim 


~ 1 arcmin 


Kessler et al. ( 1996) 


COBE/DIKBE 


100, 140 & 240 ^J,m 


~ 0.7° 


Hauser et al. (1998) 


WMAP 


23, 33, 41, 61, & 94 GHz 


53 - 13 arcmin 


Hinshaw et al. (2009) 


PMN 


4.85 GHz 


~ 5 arcmin 


Griffith & Wright (1993) 


Stoclcert 


2.72 GHz 


~ 20 arcmin 


Reif et al. (1984) 


HartRAO 


2.326 GHz 


~ 20 arcmin 


Jonas ct al^ (1998) 


Stockcrt 


1.4 GHz 


~ 35 arcmin 


Reich & Rci^h (1986) 


Parkes 


0.408 GHz 


51 arcmin 


Haslam et al. (1981) 



^ : is the angular resolution FWHM 
Spitzer program 40154 



/m(x) = Y^^{wp{x)Ip{x)/Ap{x))/Y^^Wp{x), where Wp{:x.) = 
Ap/ap, with p the label for the pth pointing and Ap{x) is 
the primary beam response. 

We also used an alternative method to restore the 
visibilities in order to check the CLEAN reconstruction. 
We applied the Vor onoi image reconstruction (VIR) from 
ICabrera et al.l (|2008| ). that is well suited to noisy data sets. 
This novel technique uses a Voronoi tessellation instead of 
the usual grid, and has the advantage that it is possible 
to use a smaller number of free parameters during the re- 
construction. Moreover, this technique provides the optimal 
image in a Bayesian sense, so the final image is unique. 

Fig. [3] shows the VIR image for LDN 1780. The point 
source was subtracted from the visibilities before the re- 
construction. The result is visibly better because the image 
negatives are less pronounced than CLEAN (Fig. [2]). The 
largest negative intensity with CLEAN is —0.007 MJy/sr 
whereas with VIR is —0.005 MJy/sr. The dynamic range is 
larger with the VIR reconstruction: 14.8 versus 11.5 from 
CLEAN. The disadvantage of the VIR reconstruction is the 
large amount of CPU time required. 

Although the VIR reconstruction is better, both tech- 
niques gives very similar results. Because of this, we trust 
the CLEAN reconstruction and used those images for the 
rest of the investigation. 



Table 3. Flux densities extracted within a 1° circular aperture 
for C Oph and LDN 1780. 



2.2 Auxiliary data 

Table |2] lists the auxiliary data used here. We used images 
form the Southern H Alpha Sky Survey Atlas (SHASSA), 
the Spitzer, IRAS, ISO, COBE and WMAP satellites, and 
low frequency data from the Parkes-MIT-NRAO (PMN) 
survey, the Rhodes/HartRAO survey, the Stockert radio- 
telescope, and the Haslam 408 MHz survey. Figures |4] and [5] 
show large fields of 3° around the two clouds from some of 
the aforementioned data. 



Freq. 
(GHz) 


Telescope / 
Survey 


91 

( arcmin) 


Flux Density 
C Oph (Jy) 


Flux Density 
LDN 1780 (Jy) 


0.408 


Parkes 


51 


55.4 ± 5.5 


39.6 ± 4.0 


1.4 


Stockert 


40 


66.9 ± 6.7 


56.9 ± 5.7 


2.3 


HartRAO 


20 


19.2 ± 1.9 


9.5 ± 1.0 


2.7 


Stockert 


20.4 


23.7 ± 2.4 


14.9 ± 1.4 


23 


WMAP 


52.8 


7.6 ± 0.1 


1.5 ± 0.1 


33 


WMAP 


36.6 


7.0 ± 0.2 


0.9 ± 0.2 


41 


WMAP 


30.6 


7.0 ± 0.2 


0.6 ± 0.2 


61 


WMAP 


21 


6.2 ± 0.4 


< 1.0 (3cr) 


94 


WMAP 


13.2 


7.2 ± 0.6 


< 2.1 (3o-) 


1249 


DIRBE 


42 


4940 ± 580 


1600 ± 90 


2141 


DIRBE 


42 


7570 ± 300 


4150 ± 320 


2997 


DIRBE 


42 


5600 ± 200 


2900 ± 160 



9 IS the angular resolution FWHM 



3 ANALYSIS 

3.1 Expected radio emission from Ha emission. 

The radio free-free emission must be accurately known in or- 
der to quantify the contribution of any dust-related excess 
emission at GHz- frequencies. One way to estimate its contri- 
bution is with the Hq emission intensity, provided that this 
line is the result of in-situ recombination and not scattering 
by dust. 

We used the continuum-corrected Ha image from the 
SHASSA survey l|Gaustad et al.l |2001^ . On the region of 
C, Oph, the image was saturated in the position of the star 
C, Oph so re-processing was necessary to get a smooth im- 
age. We masked the saturated pixels and interpolated their 
values using the adjacent pixels. Then, a median filter was 
applied to remove the field stars. The processed image was 
corre cted for d ust absorption using the E{B — V) template 
from Schlegel et al. (199^) and the extinction curve given 
bv ICardeUi et al.. (,198a ). The extinction at Ha is A(Ha) 
= 0.82 ^(V) and using Rv = A{V)/E{B - V) ^ 3.1 we 
have that A(Ha) = 2.54 E{B - V). Finally, we generated a 
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Figure 5. Maps from the auxiliary data around the extraction aperture for LDN 1780, 3° per side, in J2000 equatorial coordinates. The 
1° extraction aperture is shown. 



free-free brightness temperature map at 31 GHz using the 
relationship between Ha intensity a nd free-free brightness 
temperature presented in iPickinson et all (120031 ) assuming 
a typical electron temperature Te = 7000 K, appropriate for 
the solar neighbourhood. 



3.2 Spectral Energy Distribution fit 

In this analysis we did not use the CBI data because of the 
large flux loss on angular scales larger than a few times the 
synthesised beam. We performed simulations to estimate the 
flux loss after the reconstruction and we could only recover 
~20% of the flux. This would imply that most of the emis- 
sion from these two clouds is diffuse with respect to the CBI 
and CBI2 beam. 

The images were smoothed to a common resolution 
of 53 arcmin, the lowest resolution of the data used 
[WMAP K). We integrated the auxiliary data in a circu- 
lar aperture 1° in diameter around the central coordinates 
of the 31 GHz data. Background emission was subtracted in- 
tegrating in an adjacent region close to the aperture. This is 
necessary in the case of the low frequency radio data because 
these surveys have large b aseline uncertaintie s (see for ex- 
ampl e iReich fc Reichll2009l , or the discussion in lDavies et al.l 
Il996l ) In Table |3] we present flux densities for Oph and 
LDN 1780 in this aperture. 



a synchrotron spec tral-index {Tt oc u'^^) with /3s = —2.7 
l|Davies et al.| [l996!). The free- free power-law is fixed to the 
predicted brightness temperature at 31 GHz obtained from 
the Ha image. The s pectral index used is /3// = —0.12 
(|Dickinson et al.ll2003l 'l. The dominant source of error in the 
determination of the free-free emission through the Ha line 
is the correction for dust absorption. In Fig. [S] the dotted 
lines denote limits for the free-free law: the lower is the emis- 
sion expected from the Ha image without dust absorption 
correction and the upper dotted line is the emission expected 
with the correction (we assume that all the absorption oc- 
curs as a foreground to the Ha emission). A modified black 
body with fixed emissivity index /3 = 1.6 and 24 K fits the 
DIRBE 240, 140 and 100 fj.m dat a points. The spinning dust 
component was fitted using the lAli-Hai'moud et al.l ( 20091 ) 
models. The main parameters used in the spinning dust 
model were n{H) = 100 cm~^ and a temperature of 30 K. 
The emissivities are given in terms of the H column den- 
sity and we used Nh = 1.4 x lO^^cm"^ l|Mortonlll975l ). The 
parameters of the dust properties are reasonab le for these 
enviro nments, and were taken from iWeingartner fc Draind 
(|200ll ). 

On 1° spatial scales, the free-free emission from the Hll 
region dominates the spectrum in the cm range. The angular 
resolution of the auxiliary radio data is lower than that of 
the CBI, so we cannot make a fiux comparison on the CBI 
angular scales. 



3.2.1 C Oph 

Fig. [5] shows the Oph SED within a 1° aperture, as tabu- 
lated in Table [3] The model is the sum of synchrotron, free- 
free, thermal dust and spinning dust emissions. We used 



3.2.2 Surface brightness spectral comparison 

The PMN image (Fig. [7} has an angular resolution similar 
to that of the CBI. However this survey has been high-pass 
filtered; extended emission on scales larger than ~ 20 ar- 
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Figure 6. SED of ( Oph on 1°. The dotted red line shows the 
maximum and minimum free- free contribution expected given the 
H a emission. In green is the spinning dust emission from a model 
of lAli-Haimoud et al.l ||2009|') . 



Table 4. Surface brightness values in the peak of the 31 GHz 
image. 



Frequency 


Telescope 


Surface brightness 


Beam size 


(GHz) 




mjy/beam 


(arcmin) 


4.85 


Parkes 64m 


11 ± 4 


5.1 


31 


CBI 


25 ± 4 


7.7 



cmin is removed (| Griffith fc Wright! [l993 ). To make a com- 
parison on CBI angular scales, we simulated visibilities of 
the PMN image using the same (u, v) coverage of our CBI 
data. Since large scale modes have been removed from the 
PMN data, we exclude visibilities at {u, v) radii < 137 wave- 
lengths (corresponding to angular scales < 25 arcmin). Ide- 
ally we should consider only visibilities that correspond to 
angular scales < 15 arcmin to make a conservative analy- 
sis, but doing this excludes most of the visibilities result- 
ing in no signal. Table [4] lists the surface brightness val- 
ues at the peak of the 31 GHz image. The 4.85/31 GHz 
spectral index is at^gs — 0.43 ± 0.20. We see that al^gs is 
significantly different (at the 2.75cr level) from —0.12, the 
spectral index if the emission at 31 GHz were produced by 
optically thin free-free emission with Tg ~ 7000 K. This 
corresponds to a difference between the 31 GHz intensity 
and the value expected from a free-free power-law fixed 
to the PMN point, the non free- free specific intensity, of 
= - /4.85(4.85/31)-°-i2 = (12 ± 5) mjy beam^V 
The significance of the excess is not very high, but we note 
that this is a comparison at the peak of the 31 GHz map, 
which is coincident with the bulk of the free-free emission. 



3.2.3 LDN 1780 

At the location of LDN 1780, the Ha emission present in the 
SHASSA image is probably s cattered lig ht from the Galac- 
tic plane l|Mattila et al.ll2007l '). Recentlv. Iwitt et all (|201G| ') 
confirm this result and show that the cloud is embedded in a 
weaker diffuse Ha background, of w hich approximately half 
is due to scattered light. In contrast. Idel Burgo fc Cambresvl 
l|2006l ) state that the Ha emission is from the cloud itself and 
suggest a very high rate of cosmic rays to explain the hydro- 
gen ionisation. Whichever is the case, we can set an upper 
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Figure 7. Simulated 4.85 GHz image of Oph from PMN survey. 
Contours as in Fig. [2] 



limit to the free-free contribution using the SHASSA im- 
age. The 100, 140 and 240 fim DIRBE and WMAP 94 GHz 
points were fitted using a modified black body with fixed 
emissivity index 13 — 2. The derived temperature is 17 K. 
For the spinning dust component we used n(H) = 500 cm~"^ 
and T = 20 K. We estim ated the H col umn density using 
the extinction map from ISchlegel et al. |_|l998). This esti- 
mate used the relation from Bohlin et al.l (|1978i ') valid for 
diffuse clouds: N(}i+R2)/E{B - F)=5.8x lO^^cm'^mag-^ 
Wc found Nh = 3.5 x lO^^cm'^ 

Fig. [S] shows the fit. The pair of dotted lines sets limits 
on the free- free contribution from the Ha data. On the figure 
are also plotted 3 a upper limits to the contribution at 61 
and 94 GHz from the WMAP V and W bands. This SED 
shows an excess over optically thin free-free emission at cm- 
wavelengths. The spinning dust model is consistent with the 
SED in a 1° aperture. 



3.3 Morphological analysis 

If dust is responsible for the 31 GHz emission in these clouds, 
we expect a morphological correspondence with IR emission. 
A discussion of t he inf rared emission from dust can be found 
m iDraine fc Lil l|2007l ) and references therein. The 100 
emission is due to grains bigger than 0.01 fim that are in 
equilibrium with the interstellar radiation field at a tem- 
perature ~10-20 K. On the other hand, the mid-IR emis- 
sion traces VSGs at ~100 K. They are too hot to be in 
equilibrium with the environment so these grains are heated 
stochastically by starlight photons and, given the very small 
heat capacity of a VSG, a single UV photon increases the 
particle temperature enough to emit at A < 60 /im. 

In the SED of ( Oph the dominant contribution at 
31 GHz is free-free emission. However, inspection of the sky- 
plane images in Fig. [9] suggests that there is no correspon- 
dence between the free-free templates (Fig. [9] a,b) and the 
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Table 5. Sky correlation parameters for C, Oph. r is the linear cor- 
relation coefficient and a is the proportionality factor between the 
31 GHz image and various templates in units of (MJy/sr)~^. 
The errors are given by dispersion in the Monte Carlo simulations. 



ff 


12 fim 


25 fj,m 


60 fiin 


100 /im 


r 0.7±0.1 


0.5±0.1 


0.4±0.1 


0.5±0.1 


0.4±0.1 


a 460±100 


2.4±0.5 


2.7±0.1 


0.2±0.1 


0.4±0.1 



31 GHz contours. The CBI data seem to match better with a 
combination of free-free and IR emission. We note, however, 
that the south-eastern arm of the ( Oph cloud is slightly 
offset by ~ 3 arcmin from its IR counterpart. This offset 
is larger than the rms pointing accuracy of CBI, of order 
0.5 arcmin. 

In LDN 1780 there are clear differences among the IR 
images (Fig. I10|) . By quantifying these differences we can 
investigate which kind of dust grains (if big grains or VSGs) 
are responsible for the 31 GHz emission in this cloud. 

3.3.1 LDN 1780 

3.3.2 Sky plane correlations 

Here we investigate sky-plane cross-correlations. For this 
we used MockCBI, a program which calculates the visibil- 
ities V{u,v) of an input image I,j{x,y) given a reference 
(u, v) dataset. We computed the visibilities for the IRIS and 
free-free templates as if they were observed by the CBI with 
the (m, v) sampling of our data. We reconstructed these vis- 
ibilities in the same way that we did with the CBI data. 
Finally, we computed the correlation between the CBI, the 
IRIS and free-free templates within a square box, 30 arcmin 
per side, centred at the phase-centre of the 31 GHz data. 

To estimate the error bars in the correlations, we per- 
formed a statistical analysis. We added Gaussian noise to 
the observed visibilities and reconstructed them 1000 times. 
Then, we correlated the comparison templates with the 
noisy mock dataset. 

The histogram in Fig. Illb shows the distribution of the 
Pearson correlation coefficient of the simulated data. The 
width of the distribution gives us an estimate of the uncer- 
tainty in the correlation coefficients and also in the slope of 
the linear relation between the IR and 31 GHz data. Fig. 
Illb shows the distribution of the proportionality factor be- 
tween the emission at 31 GHz and the emission at 100 /^m. 
Table [S] lists the derived parameters. The errors are from 
the rms dispersion of the Monte Carlo simulations; they are 
conservative because of the injection of noise to the data. 

In the case of ^ Oph, there are no significant differences 
between the different IR data and the best correlation is 
with the free-free template. As we could see in the SED of 
this cloud (Fig. [5]), most of the radio emission is free-free, 
so it is not a surprise that the best correlation is with the 
free-free template. 

In LDN 1780, the correlations between 31 GHz and the 
IR data show more interesting results. Here, all the IR tem- 
plates correlate better with the CBI data than in C, Oph. The 
best match is with IRAS 60 ^m, as can be inferred quali- 
tatively from Fig. 1101 the Monte Carlo simulations confirm 
this result (Fig. 112b ). The emission at 60 fim in the diffuse 




Frequency (GHz) 

Figure 8. SED of LDN 1780. As in Fig. [6] the red dotted lines 
show upper and lower limits for the free-free contribution. The 
symbols with arrows are up per limits (3 a). In g r een is a spinning 
dust model from a model of lAli-Hai'moud et al.l ||2009| ). 

Table 6. Correlation parameters for LDN 1780. r is the linear cor- 
relation coefficient and a is the proportionality factor between the 
31 GHz image and various templates in units of (MJy/sr)~^. 
The errors are given by the dispersion in the Monte Carlo simu- 
lations. 



8 fim 


12 ^m 


25 fim 


60 fim 


100 Aim 


r 0.6±0.1 


0.5±0.1 


0.7±0.05 


0.8±0.1 


0.6±0.1 


a 5.3±1.0 


5.2±1.4 


3.7±0.9 


0.9±0.2 


0.2±0.1 



ISM is mai nly from VSGs an d a 30-40 % contribution from 
big grains (|Desert et al.|[l990l ) so our results favour the idea 
of a VSGs origin for the cm-wave radiation. Table|S]lists the 
results. 

PAHs emission at 8 & 12 fim have similar correlation 
coefficients to that of the 100 fim emission. However, the 
emission of VSGs depend on the strength of the IRF. The 
differences in morphology that appear in this cloud depend 
both on the distribution of grains within the cloud and in 
the manner that the cloud is illuminated by the IRF. In 
the next section we investigate how the IRF illuminates this 
cloud. 

3.3.3 IRF on LDN 1780 

The radio e mission from spinning VSGs i s fairly independent 
of the IRF (iDraine fc LazarianI [l998bl : lAli -Hai'moud et al.l 
l2009l : lYsard fc Verstraetd l2o'lol ). On the other hand, their 
near-IR (NIR) emission is due to stochastic heating by in- 
terstellar photons so it is proportional to the intensity of 
the radiation field and to the amount of VSGs. Therefore, 
NIR templates corrected for the IRF would trace better the 
distribution of VSGs. 

T he intensity of the radiation field Go can be estimated 
(as in lYsard et al.l I2OI0I ) from the temperature of the big 
grains Tbg in the cloud. We constructed a temperature map 
of those grains fitting a modified black body to the ISO 100 
& 200 fim images pixel-by-pixel (at the same resolution). 
With this temperature map, we calculated Go given that: 



with ^ = 2. We divided the NIR templates (8 & 12 fim) 
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Figure 9. Comparison of different templates of ^ Oph on tlie sky plane. From the left to right we show the free-free template obtained 
from the Ha image, the Stockert 2.7 GHz continuum map, the IRAS 100 fim map (which traces large dust grains), and the IRAS 12 ^m 
map (which traces very small grains). The contours show the CLEANED image of Oph. Contours are as in Fig. [2] The south-eastern 
"arm" of the 31 GHz contours has no counterpart in the free-free template; however it resembles the dust emission. 
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Figure 10. Comparison of the restored CBI2 image of LDN 1780 with IRAS templates. The radio point source NVSS 153909-065843 
has been subtracted from the CBI2 data. Contours are as in Fig. [2] Note the morphological differences between the IRAS bands. The 
arrows in the corner are perpendicular to the Galactic plane and point towards the north Galactic pole. 
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12 jam 




Figure 11. Histograms form the Monte Carlo simulations. Left: histogram of the distribution of the correlation coefficient r for the 
different templates on f Oph. The correlation with the free-free template is significantly better than the correlation with any dust 
template. Right: histogram of the distribution of a, the proportionality factor between the 31 GHz and the 100 /im images. The rms 
dispersion of these simulations are used as error bars for our results. 
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Figure 12. Histograms form the Monte Carlo simulations for LDN 1780. Left: distribution of the correlation coefficient r for the different 
IR templates. Right: histogram of the distribution of a, the proportionality factor between the 31 GHz and the 100 lira images. 



LDN 1780 



LDN 1780 




1(31 GHz) [Jy/beam] 

(a) 



d 



3, 

LO 

cv 



0.1 



-0.1 



-0.05 







0.05 



0.1 



1(31 GHz) [Jy/beam] 

(b) 



Figure 13. Pixel-by-pixel correlation of CBt data with IRAS 60 /xm and the corrected IRAS 25 /xm templates on LDN 1780. IRAS 60 /xm 
correlates better in the more diffuse regions while IRAC 25 /xm does better in the peak of the CBI image. 



by our Go map and then repeated the correlations with the 
31 GHz data. Again, to estimate error bars, we used the 
Monte Garlo simulations. The correlations in this case are 
tighter than with the uncorrected templates; here r = 0.69± 
0.04, a 2 (7 improvement. 



4 COMPARISON OF 31 GHZ EMISSIVITIES 

One motivation for this work is that the physical conditions 
in translucent clouds approach that of the cirrus clouds. 



Here, we compare the radio emission of different clouds in 
terms of their N{H). To avoid differences in beam sizes and 
frequencies observed, we choose to compare only sources ob- 
served by the CBI, at 31 GHz and scales of 4-8 arcmin. 
We also calculate an averaged column density for the cirru s 
clouds using the extinction map from ISchlegel et al] (|l998l ) 
in the positions observed bv lLeitch et aL ( 19971 ) at freguencv 
of 32 GHz and a beam size ~7 arcmin. This beam and fre- 
quency are similar to those of the CBI. 

In Table [71 we list the radio intensity at the peak of the 
CBI images, alongside a value for the column density at the 
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Table 7. Emission parameters for different clouds, all of them 
but the cirrus were observed by the CBI at 31 GHz. The second 
column is the column density, the third is the peak intensity at 
31 GHz and the fourth is the emissivity at 31 GHz. References ar e 
the following: (1) Lcitch ct al. (19 91), (2)ICasassus et al.l l l200d) . 
(Sl FCasassus et al.i (|2008|) . (4') .Castellanos et al.l 1I2OIOI') 





N{H) 


/31 ' 




Cirrus^ 


0.15 ±0.07 


6.9 ± 1 


4.6 ±2.0 


COph 


0.22 ±0.02 


9± 1 


4.1 ± 0.6 


LDN 1780 


0.45 ± 0.04 


16 ± 1 


3.5 ±0.4 


LDN 1622^ 


1.5 ±0.15 


30 ±2 


2.0 ±0.2 


p Oph3 


5.0 ±0.50 


180 ± 20 


3.2 ±0.5 


M78* 


22.8 ±0.23 


210 ± 30 


0.9 ±0.1 



Xl022 cm-2 
xlO-» MJy sr-i 
Xl0~2* MJy sr-i cm^ 
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Figure 14. Emissivity at 31 GHz versus column density for dif- 
ferent clouds observed with the CBI. The point corresponding 
to cirrus clouds wa s obtained as an average from the results of 
iLeitch et al. I lll997t) . The dashed line is the best linear fit to the 
data. 

same position and the ratio between these two quantities. 
In Fig. [14] we plot these quantities. The linear fit shown has 
a slope of 0.54 ± 0.10 MJy sr^^ cm^'' and the correlation 
coefficient is r = —0.87. There seems to be a trend in the 
direction of diminishing emissivity with increasing column 
density. Despite the large variations in column density (^^2- 
3 orders of magnitude), it is interesting that the emissivities 
of the clouds lie in a small range of ~1 order of magnitude. 
However, it is worth noting that tf the inverse relation we 
see is real, it will indicate that the anomalous emission is 
not associated with large dust grains, since their number 
increases with density, because of dus t growth. A similar 
result was obtained bv iLagachd l|2003l ). who used WMAP 
data combined with IR templates and gas tracers in the 
whole sky on angular scales of 7 deg. 



5 CONCLUSIONS 

We have presented 31 GHz data of two translucent clouds, 
C Oph and LDN 1780 with the aim of characterising their 
radio emissivities. We found an anomalous emission excess 



in both clouds at 31 GHz on angular scales of ~7 arcmin in 
C Oph and --5 arcmin in LDN 1780. 

The SED of ^ Oph on large (1°) angular scales is domi- 
nated by free-free emission from the associated Hll. Because 
of this, it is difficult to quantify the contribution from dust 
to the 31 GHz data. However, when comparing with the 
optically thin free-free extrapolated from the 5 GHz PMN 
image, we find a 2.4 a excess at 31 GHz on spatial scales of 
7 arcmin in surface brightness. 

In the SED of LDN 1780 we see an excess on 1° angular 
scales. The free-free contribution in this cloud is expected to 
be very low; the Ha emission may be scattered light from the 
IRF, so it would not have a radio counterpart. A spinning 
dust component can explain the anomalous emission excess 
in the SED. Correlations between the cm- wave data and IR- 
templates shows a trend favouring IRAS 25 & 60 ^m. The 
best match in this case is with IRAS 60 ^im although the 
peak of the CBI image is best matched by IRAC 8 /im. We 
corrected the IRAC 8 pm and IRAS 12 ^m by the IRF and 
found a tighter correlation with these corrected templates. 
In the spinning dust models, the VSGs dominate the radio 
emission. Our results support this mechanism as the origin 
for the anomalous emission in this cloud. 

The 31 GHz emissivities found in both clouds are sim- 
ilar and have intermediate values between the cirrus clouds 
and the dark clouds. The emissivity variations with column 
density are small, although we find an indication of an in- 
verse relationship, which would further support a VSGs ori- 
gin for the cm-emission. The anomalous foreground which 
contaminates CMB data comes from cirrus clouds at typi- 
cally high Galactic latitudes. Here, we see that there is not 
a large difference between the radio emissivity of cirrus and 
translucent clouds on 7 arcmin angular scales. Because of 
this similarity, translucent clouds are good places to inves- 
tigate the anomalous CMB foreground. 
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